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ABSTRACT 

The effects of merging histories of proto-objects on the angular momentum distribu- 
tions of the present-time dark matter haloes are analysed. An analytical approach 
to the analysis of the angular momentum distributions assumes that the haloes are 
initially homogeneous ellipsoids and that the growth of the angular momentum of the 
haloes halts at their maximum expansion time. However, the maximum expansion time 
cannot be determined uniquely, because in the hierarchical clustering scenario each 
progenitor, or subunit, of the halo has its own maximum expansion time. Therefore 
the merging history of the halo may be important in estimating its angular momen- 
tum. Using the merger tree model by Rodrigues & Thomas, which takes into account 
the spatial correlations of the density fluctuations, we have investigated the effects of 
the merging histories on the angular momentum distributions of dark matter haloes. 
It was found that the merger effects, that is, the effects of the inhomogeneity of the 
maximum expansion times of the progenitors which finally merge together into a halo, 
do not strongly affect the final angular momentum distributions, so that the homoge- 
neous ellipsoid approximation happens to be good for the estimation of the angular 
momentum distribution of dark matter haloes. This is because the effect of the differ- 
ent directions of the angular momenta of the progenitors cancels out the effect of the 
inhomogeneity of the maximum expansion times of the progenitors. 

The contribution of the orbital angular momentum to the total angular momentum 
when two or more pre-existing haloes merge together, was also investigated. It is shown 
that this contribution is more important than that of the angular momentum of diffuse 
accreting matter to the total angular momentum, especially when the mergers occur 
many times. 
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1 INTRODUCTION 

Galaxy formation is one of the most important problems in 
astrophysics. In general, it includes different physical pro- 
cesses such as the evolution of cosmological fluctuations, 
star formation, heating processes of baryonic gas from super- 
novae, dynamical and chemical evolution of gas, etc. Hence 
it is difficult to attack the problem of the formation and 
evolution of galaxies in a way that connects all of the com- 
plicated physical processes. 

One of the important keys to understanding galaxy for- 
mation is the angular momentum of dark matter haloes. For 
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example, there is a strong correlation between the angular 
momentum and the morphology of galaxies. The mean value 
of the angular momentum of ellipticals is much smaller than 
that of spirals. Thus we may expect to understand the origin 
of the morphological distinction from analysis of the angu- 
lar momentum. In the standard model, it is considered that 
dark matter dominates in our Universe and that galaxies 
and clusters of galaxies have been formed by the gravita- 
tional growth of small initial density fluctuations. The dark 
matter has collapsed and virialized by self-gravitational in- 
stability into objects which are called 'dark matter haloes' 
or 'dark haloes'. The large haloes are generally considered to 
have been formed hierarchically by the clustering of smaller 
ones (so-called 'hierarchical clustering'). During this pro- 
cess, each dark halo obtains its angular momentum by tidal 
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forces (Hoyle 1949; Peebles 1969; White 1984; Barnes & Efs- 
tathiou 1987). Also, some properties of the baryonic compo- 
nent of galaxies e.g. morphologies, surface brightness and so 
on, have been discussed as a function of angular momentum 
of dark haloes (e.g., Dalcanton, Spergel & Summers 1997; 
Jimenez et al. 1997; Mo, Mao & White 1998). 

Doroshkevich (1970) and White (1984) found that the 
angular momentum evolves proportionally to t (t is the 
cosmic time) during the linear regime, considering the ini- 
tially homogeneous density ellipsoid as a proto-object in 
the Einstein-de Sitter universe. This prediction has been 
confirmed by the N-body simulations by White (1984) and 
Barnes & Efstathiou (1987). It should be noted that White 
(1984) did not consider whether the dark matter in the re- 
gion under consideration would in fact grow into a collapsed 
object or not in his analytical prediction. 

In the analysis of the evolution of the angular momen- 
tum of a proto-object, it is assumed that the dark mat- 
ter around a density peak maximum collapses into a halo 
(Hoffman 1986, 1988; Ryden & Gunn 1987; Ryden 1988; 
Caimmi 1989; Quinn & Binney 1992 and Eisenstein & Loeb 
1995a, b). The angular momentum distributions were de- 
rived analytically by Heavens & Peacock (1988), and Cate- 
lan & Theuns (1996a, b). They used the peak formalism of 
Gaussian random fields (Peacock & Heavens 1985; Bardeen 
et al. 1986), in order to analyse statistically a density pro- 
file around the density peak maximum. Catelan & Theuns 
also considered the case of non-Gaussian initial conditions 
(Catelan & Theuns 1997). Susa, Sasaki & Tanaka inves- 
tigated the angular momentum distribution by using the 
Press-Schechter formalism (Press & Schechter 1974) instead 
of the peak formalism. 

It should be noted that when the distribution of the 
angular momentum was derived, there was assumed that an 
'object' is the initially homogeneous density ellipsoid and 
that the acquisition of the angular momentum halts at the 
maximum expansion time of the object (Heavens & Peacock 
1988; Catelan & Theuns 1996a, b). The maximum expansion 
time is usually estimated from the averaged density contrast 
of the object by assuming spherically symmetric collapse. 
In spherical collapse, the averaged linear density contrast of 
the object reaches to 1.05 at the maximum expansion time 
(e.g. Peebles 1993). However, in the hierarchical clustering 
scenario, small objects collapse firstly. Some of the small ob- 
jects may then merge together, leading to the present-time 
large halo. It is general in this scenario that the acquisi- 
tion of the angular momentum of a dominant proto-object 
of a halo may halts earlier than the maximum expansion 
time estimated through the averaged density contrast of the 
present-time halo. Hence, if we believe in the merging his- 
tory of the proto-objects, then we should take into account 
the halting time of each proto-object, a contribution of the 
orbital angular momentum of the merging proto-object, and 
the angular momentum of the accreting matter. 

Merging history models of dark haloes have been used 
for the semi-analytical methods of galaxy formations (Cole 
& Kaiser 1988; Kauffmann & White 1993; Rodrigues & 
Thomas 1996; Roukema et al. 1997). Cole & Kaiser devel- 
oped the block model, considering the one-point distribu- 
tion function of density fluctuations. Therefore spatial cor- 
relations of the density fields were not taken into account. 
Kauffmann & White (1993) extended the Press-Schechter 



formalism and constructed the merging histories by coupling 
this formalism with the Monte Carlo method. The extension 
in question is mainly based on the results of Bower (1991). 
Kauffmann & White (1993) have utilized the one-point dis- 
tribution function as well. On the other hand, Yano, Na- 
gashima & Gouda (1996) have shown that the spatial cor- 
relation of the density field is important when calculating 
mass function of dark haloes. The analytically derived mass 
function accounting explicitly for the two-point correlation 
function, differs significantly from that found by the Press- 
Schechter mass function. So we believe that the merging 
histories should also include spatial information. Rodrigues 
& Thomas (1996) constructed the merging history model 
(which we call 'the merging cell model', for convenience), 
including the information on the spatial correlations. In the 
merging cell model, the random Gaussian density fluctua- 
tion field is realized on spatial grids as in the construction of 
the initial conditions of N-body simulations. Therefore this 
model includes, naturally, information on the spatial corre- 
lation. Then, by finding the collapsed region whose linear 
density contrast, 5, averaged over the region is 8 C = 1.69 (8 C 
is the critical density contrast for collapse, see, e.g., Peebles 
1993), we can construct merger trees. We expect that this 
model is more realistic and useful in describing the galaxy 
formations although a spherical collapse of each block is 
assumed. Nagashima & Gouda (1997) have analysed this 
merging cell model by comparing the mass functions given 
by the merging cell model with the mass functions by Yano, 
Nagashima & Gouda (1996). They found that both types of 
the mass functions are in good agreement. 

In this paper, using the merging cell model, we calcu- 
late both the density fields and the velocity fields in order to 
estimate the angular momentum distribution of dark haloes. 
In a similar way to the previous analytical work (Heavens 
& Peacock 1988; Catelan & Theuns 1996a, b), the initial 
angular momentum of each halo is calculated through the 
velocity field, then the angular momentum of each halo is 
evolved proportionally to the cosmic time t. Since the over- 
density of each halo is known, we can estimate its maximum 
expansion time. By using this time, we compute the an- 
gular momentum of the halo. When a merger between the 
haloes occurs, we investigate the spins and the orbital an- 
gular momentum of these haloes around the centre-of-mass 
of a larger new common halo, and the angular momentum 
of the matter being accreted into the new halo. The sum 
of these angular momentum components is the spin of the 
new halo. This process is repeated when the mergers occur. 
It should be noted that the block model and the extension 
of the Press-Schechter formalism cannot deal with angular 
momentum, because these models ignore the velocity field 
structures. 

To compare our method with the previous work, we 
also calculate the angular momentum of homogeneous den- 
sity haloes by smoothing density contrasts in each region 
occupied by haloes at the present epoch. In this way we in- 
vestigate the effects of the merging histories by showing the 
angular momentum distribution of dark haloes. 

Our aim is to clarify the physical effects of the merging 
histories of dark haloes, namely, the effects of the differ- 
ence in the maximum expansion times and the difference in 
the directions of the angular momenta of the proto-objects, 
that is, the subunits of the haloes. An advantage of our 
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method is that it makes available a qualitative analysis of 
the acquisition of angular momentum. Of course, we must 
perform N-body simulations in order to understand the pro- 
cess quantitatively (Efstathiou & Jones 1979; Barnes & Ef- 
stathiou 1987; Warren et al. 1992; Lemson & Kauffmann 
1997), but our calculations require less computing time than 
do N-body simulations. Therefore, we can investigate many 
model parameters for a short time. We believe that the semi- 
analytical methods are complementary to N-body simula- 
tions. 

In Section 2, previous analytical work is reviewed 
briefly. In Section 3, the merging cell model is character- 
ized. In Section 4, the method of calculating the angular 
momenta of haloes in the merging cell model is shown. In 
Section 5, we consider the angular momentum distribution. 
We show the merger effects as well as the role of the orbital 
angular momentum of the merging objects in their common 
halo. Section 6 is devoted to conclusions and discussion. 



2 TIME DEPENDENCE OF ANGULAR 
MOMENTUM 

In this section, we review the evolution of angular momen- 
tum according to the White's method (White 1984). The no- 
tations are almost the same as in Catelan & Theuns (1996a). 

The angular momentum L of an object occupying the 
volume V in the Eulerian space is defined as 

L(t)= f p(r)(r-r )x(v-v )dr, (1) 
Jv 

where p is the density, r is the Eulerian coordinate, v is the 
velocity, and the subscript denotes the centre-of-mass of 
the volume V. This transforms to the Lagrangian descrip- 
tion, 

L = y?7oa 2 [q-qc) + DV (<t> - <t>o)] xDV (<f) ~ <f> )dq 

= J 77 a 2 I)(q-qo)xV(^-^o)dq, (2) 

with r/o = pa 3 , where 770 is the comoving mean density of 
the universe, p is the mean density of the universe, a is 
the cosmic scale factor, q is the Lagrangian coordinate, <j) is 
proportional to the gravitational potential in the Einstein- 
de Sitter universe, D is the linear growth factor (D — a 
in the Einstein-de Sitter universe), and V denotes the La- 
grangian volume corresponding to V. The dot represents 
time derivative, d/dt. In the above transformation, we use 
the Zel'dovich approximation (Zel'dovich 1970), 

r(t,q)=o(t)[q + I>(t)V#q)]. (3) 

Thus we find the growth rate of the angular momentum as 

L oc aD oc t (4) 

in the Einstein-de Sitter universe. 

The growth of the angular momentum halts at the max- 
imum expansion time of the object. In the previous work, 
this maximum expansion time is estimated via the density 
contrast at a centre of the object with the Lagrangian vol- 
ume T, which is smoothed with a filtering scale correspond- 
ing to the mass of the object. Using the spherical collapse 



approximation, the maximum expansion time t m is obtained 
by t m = to(1.05/5) 3 / 2 , where to is the present time and <5 is 
the density contrast at the centre of r. 

Note that the above discussion is justified only when 
the Zel'dovich approximation is valid anywhere in F until 
tm- However, in hierarchical clustering scenarios, collapsed 
haloes cluster and grow up gradually, so that the Zel'dovich 
approximation in such pre-existing haloes is broken earlier 
than t m of the final collapsed halo with V. Therefore, we 
need to investigate the case in which the growth of angular 
momenta of the pre-existing haloes halts earlier than the 
maximum expansion time of the final collapsed halo which is 
formed via mergers of the pre-existing haloes. In this paper, 
a progenitor means such a pre-existing halo. 



3 MERGING CELL MODEL 

Let us describe the general aspects of the merging cell model 
in accord with Nagashima & Gouda (1997). 

First, the random Gaussian density field is realised in 
a periodic cubical box of side L. The Fourier mode of the 
density contrast S[= (p — p)/p] obeys the following proba- 
bility distribution function for its amplitude and phase. In 
the random Gaussian distribution, (Bardeen et al. 1986), 

P(|5 k |,0 k )d|«5 k |^ k = ^iexp|-M|d|«5 k |^i, (5) 

where <j> k is the random phase of 8 k , 8 k = \5 k \ exp (i(j> k ), and 
P(k) is the power spectrum. Then, the density contrast at 
each grid ('cell') is given by the Fourier transform, 

5(q) = ^ C s ^^ k > (6) 

where k c is the cut-off wavenumber corresponding to the cell 
size. 

Next, the fluctuations of the various smoothing levels 
are constructed by averaging the density fluctuations within 
cubical blocks of side 2, 4, . . . , L^, ox - At each smoothing 
level, displacing the smoothing grids by half a blocklength 
along the coordinate axes, eight sets of overlapping grids are 
constructed to centre approximately the density peak within 
one of them. 

Then, the density fluctuations within blocks and cells 
are combined into a single list and ordered by decreasing 
density. The fluctuations are investigated from the top of 
this list. The following criteria decide whether a block or a 
cell can collapse. Note the terminology that a halo is a block 
or a cell that has already collapsed, and an investigating 
region is a block or a cell whose linear density contrast is 
just equal to 8 C at the reference time. 

(a) If the investigating region includes no haloes, its can 
collapse and can be identified as a new halo. 

(b) When the investigating region includes a part of a 
halo, if the overlapping region has at least half of the mini- 
mum of the masses of the halo and the investigating region, 
then the investigating region can collapse. This is the crite- 
rion of collapse of the investigating region and merging for 
the overlapped haloes. We call this criterion as the overlap- 
ping criterion. 

(c) If the investigating region has half or more of its mass 
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contained in exactly one pre-existing halo and the overlap- 
ping criterion (b) is satisfied for any overlapped haloes then 
they can be merged together as part of a new structure. This 
is the criterion for linking of haloes. 

These criteria were chosen by Rodrigues & Thomas 
(1996). We introduce the overlapping parameters x and y, 
according to Nagashima & Gouda (1997). The parameter x 
means the ratio of the mass of the overlapping region to the 
lesser of the masses of the halo and the investigating region. 
Then, the Rodrigues & Thomas's criterion corresponds to 
x — 1/2. The parameter y quantifies the linking criterion. 
Using this parameter y, we change the criterion (c) as fol- 
lows: If the investigating region has y times (x < y < 1/2) or 
more of its mass contained in exactly one pre-existing halo 
and the overlapping criterion (b) is satisfied for any over- 
lapped haloes then they can be merged together as part of 
a new structure. 

In the case of x = y = 1/2, it is shown by Na- 
gashima and Gouda (1997) that the mass functions given 
by the merging cell model are in agreement with those of 
Yano, Nagashima and Gouda (1996), whose formula implic- 
itly includes an overlapping criterion which corresponds to 
x = y — 1/2 as an assumption. Thus this agreement is lost 
in the case of x = 1/8 and y — 1/2. However, we need more 
detailed investigations by, e.g., N-body simulations, in order 
to decide realistic values of x and y. Therefore we use two 
parameter sets x = y = 1/2, and x — 1/8 and y = 1/2 in 
this paper. 



4 METHOD OF CALCULATION OF 
ANGULAR MOMENTUM 

In this section we describe the method of calculation of the 
angular momentum of haloes. 

From the Fourier component <5k obtained by the Monte 
Carlo method using eq.(^|) and the Poisson equation 
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Figure 1. Example of generation and merging history of haloes. 
The first generation haloes include no progenitors. The final halo 
belongs to the fourth generation, and is formed by the merger 
of the third generation halo and the second generation halo. We 
omit the accreting matter in this figure. 



contrast of the new common halo in order to estimate the 
maximum expansion time. We estimate the relation between 
the linear density contrast at the present epoch and the max- 
imum expansion time t m using the spherical collapse model, 



tm — to 



/ 1.05 
V S 



(10) 



Now we define the generation of haloes as follows (see 
Fig.[j]) . A collapsed block which has no progenitors is defined 
as a first generation halo (k = 1). The fc-th generation halo is 
defined as a halo which includes the (Je— l)-th generation pre- 
existing haloes whose number specifying their generations 
is the largest of its progenitors. If the new halo is a fc-th 
generation halo (k > 2), these three components of angular 
momentum are expressed as 



we can estimate the velocity v and potential <j> at each cell, L sp i„ = S ' L fc ^ [1 < fc'(i) < k — 1] 



(qj = 



v 



v(q) 



(27r) 3 £>0 

. a D V 



°k ikq 
P 

k 2 



d 3 k, 



(2-kYDoJ k 2 



kSk e^d 3 k, 



(8) 



(9) 



^M;(q; - qcAf) x (vj - Vcm) 



1.05\3 



6 k 



(11) 



,(12) 



where ao and Do are the scale factor and the growing mode 
at the present epoch. 

When a halo collapses, we divide this new halo's angular 
momentum into three components, which are, (1) the sum 
of the spin angular momenta of pre-existing haloes L sp i n , 
(2) the sum of the orbital angular momenta of pre-existing 
haloes L or 5, and (3) the angular momentum of the mat- 
ter accreting into a new halo, L acc . In order to calculate 
the above quantities, we need the mass of the pre-existing 
haloes, and the positions and velocities of the center-of- 
mass of both the pre-existing haloes and the new common 
halo. The growth of the angular momenta of the orbital and 
accreting matter is proportional to t, according to eq.(Q). 
These components grow until the maximum expansion time 
of the new common halo. Therefore, we need the density 



M ceU ^(qj - qcjf) x (v 



ma* 



where the subscripts i and j stand for the isolated pre- 
existing haloes just before the new common halo col- 
lapses and cells of accreting matter respectively, so that 

k' (i) 

L spin i , Mi, qi and v» are respectively the spin, mass, po- 
sition and velocity of the centre-of-mass of the i-th progen- 
itor. The superscripts k and k'(i) show the generation of 
haloes. The isolated fc'(i)-th generation pre-existing haloes 
[k'(i) > 1] merge together into the fc-th generation halo. 
The subscript CM denotes the centre-of-mass of the new 
halo, M ce u is a mass of one cell, and 8 k is associated with 
a new collapsed fc-th generation halo, because the angular 
momenta of the orbital and accreting matter grow until the 
maximum expansion time of the new collapsed halo. Finally, 
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Figure 2. Schematic description of dark halo formation in- 
cluding two pre-existing haloes. The halo labelled by the su- 
perscript k, which is surrounded by the ellipsoid r fe , is formed 
through the merger of the two progenitors labelled by k — 1, 
which are surrounded by the ellipsoids 1 and T^ 1 , and 
the accretion of matter within T . These two progenitors have 
their spins I>J^J, 1 and Lj^"* 2 > and their orbital angular mo- 
menta around the centre-of-mass of the new common fc-th gen- 
eration halo L or (, i and L or f,,2- These orbital angular momenta 
and the accreting matter's angular momentum L acc grow until 
the maximum expansion time of the new common halo. Thus 
the angular momentum of the new common halo is obtained by 

j n — Tj S p2n T L or b -f- L acc . 



we obtain a total angular momentum L* pi7i of the new col- 
lapsed halo by summing the above quantities, 

+ L orb + L acc . (14) 

If this new collapsed fc-th generation halo will merge 
again into a collapsed halo at a later epoch, Lj pin is used 
as L* pin j in the estimation of the angular momentum L fc+1 
of the later collapsed (k + l)-th generation halo. Note that 
when a first (k — 1) generation halo collapses, this halo has 
no progenitors, so this halo is formed by only the accreting 
matter. Therefore the angular momentum of this halo is es- 
timated by only eq.(^), i.e., L* pin = L acc . This estimation 
is the same as in the homogeneous collapse, eq.(^). 

When considering a homogeneous case as in the previ- 
ous analytical work (Heavens & Peacock 1988; Catelan & 
Theuns 1996a, b), we first calculate the regions of collapsed 
haloes as described in Section 3. Then we average the den- 
sity contrasts of cells in each halo, and obtain the averaged 
density contrast 8 avg . Finally, using only eq. (|l3|) , we obtain 
the angular momentum of each halo. In this equation, 8 k 
corresponds to the above averaged density contrast 5 a v g - 

In Figj^, we show the schematic description of dark halo 
formation including two progenitors. In this case, two (fc— 1)- 
th generation progenitors merge together into a new fc-th 
generation larger halo with accreting matter. So the spin 
component of the new fc-th generation halo is obtained as 
Ls P in = Lspi^i +Lj~^ 2 , and the total angular momentum, 



namely, the spin of the new fc-th generation halo is L* pin = 



5 RESULTS 

In this section we show final distributions of the spin angular 
momentum |L| and the mass- weighted angular momentum 
L|/Af 5//3 , which are similar to those used by Heavens & Pea- 
cock (1988) and Catelan & Theuns (1996a). We assume that 
linear density fluctuations obey a random Gaussian distri- 
bution with a power-law power spectrum P(k) tx fc™, where 
n — and —2. We consider only the Einstein-de Sitter back- 
ground universe (f2 = 1, A = 0). The calculating box size is 
Lu QX = 128. The normalization of the power spectrum is 
such that the variance of the density contrasts equals unity 
at a block with eight cells. The overlapping parameters are 
set as x — y = 1/2, and x = 1/8 and y = 1/2. 

The number of haloes finally obtained is about 1.6 x 10 4 
for n — and about 9 x 10 3 for n — —2 in one calculating 
box. The following results are averaged over five realizations, 
therefore our investigations are statistically robust. 

In all figures, the displayed values of the mass, the an- 
gular momentum and the time are normalised. The time t 
is normalized by the present time to- Physical values for the 
mass and the angular momentum at the present epoch may 
be obtained through M v h ya = m,M and T-i p h y s = m*£*L, 
where m« = 1.4 x 10 14 Af e and I, = 3 x 10 6 km s _1 kpc 
for the Hubble parameter h — 0.5, which is normalised 
by 100 km s _1 Mpc _1 (m» corresponds to the mass of one 
cell and U is determined by the resolution of the calcula- 
tion). Note that these physical values are scaled as m*(z) = 
m» (1 + z }- 6 / C i + n ) because of the scale-free power spectrum. 
For example, m» ~ 3 x 10 9 Mq at z — 5 in the case of 
n = -2. 

In the following, we consider three cases in which to 
estimate the angular momentum. One is the merger case, 
in which L is given by eq. (p^|) . Second is the non-orbital 
case, in which we remove the component of orbital angu- 
lar momentum, L or (,. Third is the homogeneous case. In the 
homogeneous case, we estimate a region occupied by each 
halo, then we average the density contrasts in the region; 
after that, we estimate the maximum expansion time via 
the averaged density contrast by using eq.(|l(]). Thus each 
cell which belongs to a halo has the same maximum ex- 
pansion time, that is, the halo is treated as a homogeneous 
density object in the homogeneous case. By using this max- 
imum expansion time, we obtain the angular momentum of 
the final halo through the velocity field, as in eq.([l3]). This 
procedure is almost the same as that in the previous work 
(Heavens & Peacock 1988; Catelan & Theuns 1996a, b). 

In all cases, we consider only haloes with mass larger 
than nine cells, because of the numerical resolution. Also, 
we consider only final haloes with an averaged density con- 
trast larger than 1.69 in the homogeneous case, because each 
collapsing block must satisfy the criterion of the spherical 
collapse. 

5.1 Mass dependence of angular momentum 

First, we show the mean angular momentum against the 
mass of the haloes in the case of x = 1/2 (FigM and x = 1/8 
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Figure 3. Mass dependence of mean angular momentum in the 
case of x = 1/2. upper panel: n=0, lower panel: n=-2. The dashed 
lines show the homogeneous case (HOM). The solid lines show the 
merger case (MRG). The dotted lines show those without orbital 
angular momentum (the non-orbital case; NOR). The thick solid 
lines show the lines proportional to L tx M 5 ' 3 . 

(FigJiJ). The horizontal axis denotes the mass, that is, the 
number of cells of each halo. The vertical axis denotes angu- 
lar momentum averaged over each mass bin. The solid lines 
correspond to the merger case. The dotted lines and dashed 
lines correspond to the non-orbital case and the homoge- 
neous case, respectively. The short thick solid lines denote 
L tx M 5/3 . The minimum mass in the figures corresponds 
to a halo with eight cells (log 8 ~ 0.9), and the values of the 
mean angular momentum in the minimum mass are unphys- 
ically small because of the numerical resolution. In the ho- 
mogeneous case, there are no data for larger-mass because 
the averaged density contrasts in large haloes are smaller 
than 1.69. 

The deviation from the thick line L oc M 5 ^ 3 probably 
results from the boundary effect of objects, caused by the 
artificial blocks and cells. We will show this boundary effect 
in a future work. 

In the case of x = 1/2 (Fig.[|), the differences among 
the lines are rather small, but we find tendencies: the an- 
gular momentum in the merger case is smaller than in the 
homogeneous case, and the angular momentum decreases 
slightly when the orbital angular momentum component is 
removed. We find that orbital angular momentum does not 
play an important role in the case of x = 1/2. In the case of 
x = 1/2, haloes grow almost solely through accretion events; 
large blocks with 8 3 or more cells should collapse in order to 
merge two or more pre-existing haloes. However, the proba- 
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Figure 4. The same as Fig.j^ but for x = 1/8. upper panel: n=0, 
lower panel: n=-2. 

bility that such large blocks would collapse, including two or 
more haloes and satisfying the criterion (c), is small. Thus 
merger events between two or more pre-existing haloes oc- 
cur rarely. This is the reason why orbital angular momentum 
does not play an important role in the case of x = 1/2. 

On the other hand, when i = 1/8 we can see that the 
difference between the merger case and the non-orbital case 
is large, about one order of magnitude at log M ~ 3.5, while 
the difference between the merger case and the homogeneous 
case is not pronounced (Fig^). Under these conditions, the 
merger events occur easily. Since the orbital angular mo- 
mentum can grow until a relatively later time, the role of 
the orbital angular momentum becomes important. 

In both cases, we find that the mass dependence of the 
homogeneous case is almost the same as that of the merger 
case, while the contribution of the orbital angular momen- 
tum to the total angular momentum depends on the over- 
lapping parameter x. These properties are almost indepen- 
dent of the spectral index n. It seems to be strange that 
the merger and the homogeneous cases lead to the same re- 
sults independently of the value of x. We will investigate this 
problem by considering the angular momentum distribution 
in detail in the next subsection. 



5.2 Angular momentum distributions 

To investigate the angular momentum distribution in detail, 
let us consider the distributions of the angular momentum L 
in Fig.^a (x = 1/2) and in Fig.^a (x = 1/8), and the mass- 
weighted angular momentum L/M 5 ^ 3 in Fig.^b (x = 1/2) 
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and in Fig.|6p (x = 1/8). Each figure includes the merger 
case and the homogeneous case. The solid lines show the 
distributions for the merger case, and the dashed lines the 
homogeneous case. The thick lines correspond to the spectral 
index n = 0, and the thin lines mark n = —2. 

In the case of x = 1/2, we can see that the distributions 
for the homogeneous case are almost the same as those of 
the merger case. The degree of the difference between them 
is the same order of magnitude as the difference between the 
mean values shown in Fig.[| This result suggests that the 
final distribution of the angular momentum does not change 
very much when considering the merger effect. 

For x = 1/8, it seems that the distributions of L in the 
merger case become wider especially in the case of n = —2 
(FigJqa). It should be noted that haloes with averaged den- 
sity contrast smaller than 1.69 are removed in the homo- 
geneous case, so the difference between the homogeneous 
and the merger cases is apparently large at higher angular 
momentum. In Fig.^3, we do not need to take care of this 
removal effect in the homogeneous case, because the mean 
angular momentum is nearly proportional to M 5 ^ 3 (see the 
thick solid line in Fig.^|). The distribution for the homoge- 
neous case is also almost the same as that for the merger case 
when x = 1/8. We find that the merger effect does not affect 
the final angular momentum distribution, independently of 
x. 

Still, the difference between the merger case and the ho- 
mogeneous case is very small. This fact is considered below 
in more detail. 



5.3 Distributions of maximum expansion time 

In this subsection, first we investigate the distribution of a 
product of mass and the maximum expansion time of each 
halo, because the absolute value of angular momentum de- 
pends strongly on the maximum expansion time. This quan- 
tity is represented as 



M(tm) 



M ce ut 



(15) 



where t mt i is the first maximum expansion time of the cell 
labelled by i, which is defined below in each model. The 
sum is carried out in the region of each halo. So the angle 
brackets mean the average in the region of each halo. We 
distinguish the merger case from the homogeneous case by 
superscripts 'MRG' and 'HOM', e.g., (t^ RG ) and (t^ OM ), 
because values of t m ,% are different in the case of the merger 
and the homogeneous cases. In the homogeneous case 
of a cell is estimated from the density contrast averaged over 
the region of each halo, therefore cells which belong to the 
same final halo have the same maximum expansion time. 
However, in the merger case t™ G is estimated from the 
density contrast of a block to which the cell belongs. 

We show the distributions of M{t m ) in Figjija (x = 
1/2) and in Fig|a (x = 1/8), and (t m )/M 2/:i in FigQa 
(x — 1/2) and in Fig.^D (x = 1/8), in the merger and the 
homogeneous cases, respectively, because the dependences 
on the maximum expansion time L and L/M 5//3 in FigsJB] 
and|^are the same as M(t m ) and (t m ) /M 2 ^ 3 . The solid lines 
show the distributions for the merger case, and the dashed 
lines represent the homogeneous case. The thick lines and 
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Figure 5. (a) Angular momentum distributions in the case of 
x = 1/2. The thick lines denote the case of n=0, and the thin 
lines n=-2. The solid lines denote the merger case, and the dashed 
lines homogeneous case, (b) Mass-weighted angular momentum 
(L/M 5 / 3 ) distributions in the case of x = 1/2. Types of the lines 
are the same as (a). 



the thin lines mark the spectral index n = and n = —2, 
respectively. 

In the case of x = 1/2, the differences between the 
merger case and the homogeneous case are about 0.3 order of 
magnitude at log(M{i m )) ~ 1 — 2 (Fig.^,) and at the range 

of \og{M{tm)) ~ -2.5 1.5) (Fig.03). These differences 

are large as compared with the distribution of the angular 
momentum shown in Fig. 5| This is caused by the numerous 
cells which halt the growth of the angular momentum earlier 
in the merger case than in the homogeneous case, especially 
for large haloes. The result in the case of x — 1/8 is the same 
as for x = 1/2. Therefore, we must explain the reason why 
the distribution of the angular momentum is not changed 
much by the merger effect, as compared to the distribution 
of the maximum expansion time. 

We should consider the effect of the direction of the an- 
gular momentum. If the growth of the angular momenta of 
the proto-objects with a direction close to the direction of 
the final angular momentum halts later than the growth of 
the angular momenta of the other proto-objects with an op- 
posite direction, in a region which becomes a present-time 
halo, the total angular momentum of the halo can become 
larger in spite of the averaged maximum expansion time 
(t^( RG ) being small. Thus we investigate the maximum ex- 
pansion time, weighted by the direction (which will be de- 
fined below). 

Defining the time-independent vector part of the angu- 
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Figure 6. The same as Fig.g, but for x = 1/ 

lar momentum of a cell labelled by i as 
Li ec = M ce ii(cn - c\cm) x (vi - vom), 



(16) 



the absolute values of the angular momentum for the merger 
case and for the homogeneous case are as follows 



Ej vec.MRG 



J2i Mcellt 



HOM 
tn.i 



M 



E L " f 



(17) 



respectively. Considering the merger case, defining Of as the 
angle of the direction of the final angular momentum, we 
can transform eqs. (JlTf) and (n~8[) into 



= 22 i L r c i cos(6»i - Of)t™f 



iL^I-^lLricos^-^)^)- 



(19) 
(20) 



So we investigate the direction-weighted maximum expan- 
sion time instead of (t m ), 



(tm) 

Then 



| L MRG| „ \Lr C \ C0S(^ - 9f){t m ) c 



(21) 



(22) 



Here we have used the assumption that the absolute value 
of L" ec is statistically independent of the direction of L" ec . 



Figure 7. Distributions of product of mass M and maximum 
expansion time (t m ) in the case of a; = 1/2. (a) distributions of 
M (t m )- (b) (t m )/M 2 / 3 . The solid lines show the distributions for 
the merger case. The dashed lines show those for the homogeneous 
case. 



Thus if (t; 
IT MRG\ 



') ~ (tm)cosfl it is expected that |l/ 



We show the relation of the direction-weighted maxi- 
mum expansion time (i m )cosS, the non- weighted maximum 
expansion time (i* fflG ), and the averaged maximum expan- 
sion time (t m OM ) in Fig.| {x = 1/2) and in Fig.fu] (x = 1/8). 
Each dot corresponds to a halo. The left panels show results 
for the spectral index n — and the right panels for n = —2. 
We show the comparisons of (t^ OM ) with (i m ) CO s in the up- 
per panels, and with (tJ, fflG ) in the lower panels. The solid 
straight lines have a unity slope. 

The difference of the extents of the distributions of dots 
in the horizontal axis between the case of n = and n — —2 
is a result of the difference of the power spectrum of den- 
sity fluctuation field. The mass dependence of the maximum 
expansion time is 



3/2 
oca,'. = 



/ 1.05 
V 5 



3/2 



AT 



(23) 



thus the distribution of dots in the case of n = spread out 
wider than those in n = — 2. 

In the lower panels, there are most of dots below the 
solid lines, as we have already found in Figs.^ and [| How- 
ever, in the upper panels, the dots are distributed in a wider 
range in the vertical axis than those in the lower panels. This 
result shows that in some merging histories, the growth of 
angular momentum with the opposite direction to the final 
one happens to halt earlier than that with the same direc- 
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Figure 8. The same as Fig.]?] but for x - 



Figure 9. Distributions of maximum expansion time (t m ) in the 
case of x = 1/2. The solid lines denote those with a unity slope, 
(a) Distributions of (i^ OM ) (the homogeneous case) v.s. (t ra ) cos j 
in the case of n = 0. (b) Distributions of ((^ OM ) (the homoge- 
neous case) v.s. {t^ RG ) (the merger case) in the case of n = 0. 
(c) The same as (a), but for n = —2. (d) The same as (b), but 
for n = —2. 



tion of the final one, and so some haloes satisfy the condition 
(tm)cos9 > (t^ RG ). Then <i m )cos9 ~ (C° M ) while the dis- 
persion of (i m )cos0 is large. 

This is the main reason why the distribution of the an- 
gular momentum in the merger case is almost the same as 
that in the homogeneous case. Although (i* MG ) is smaller 
than (im OM ), proto-objects having an angular momentum 
with the same direction as the final one have a larger max- 
imum expansion time t^f G , and those with the opposite 
direction have smaller maximum expansion time. So the fi- 
nal angular momentum becomes almost the same as that 
for the homogeneous case. Therefore both distributions are 
very much alike. 



5.4 Orbital angular momentum 

In order to investigate the contribution of the orbital angu- 
lar momentum to the total angular momentum, in Fig.|ll| 
we show the case in which we remove the component of the 
orbital angular momentum (the non-orbital case) in the case 



of x = 1/2. The solid lines show the total angular momen- 
tum (Lspin + L acc + li or b) in Fig.[i"l|a and those divided by 
M 5//3 in Fig.[ll|b for the merger case, and the dashed lines 
show the total angular momentum without the orbital angu- 
lar momentum component, that is, L sp i n + L acc in Figjllja 
and those divided by M 5/3 in Fig.[n]b. The thick lines de- 
note the spectral index n = 0, and the thin lines n — —2. In 
the case of x = 1/2, reflecting the small differences in Fig.^, 
the differences between the merger case and the non-orbital 
case are very small. As already mentioned in Section 5.1, 
in the case of x — 1/2, haloes grow almost solely through 
the accretion events, and the merger events between two or 
more pre-existing haloes occur rarely. Thus the contribution 
of the orbital angular momentum is small. 

In Figjl^, we show the same results as in Figjll] but 
for the case of x = 1/8. In the upper panel, in contrast 
to Fig.prja, the differences at higher angular momentum are 
larger than those at lower angular momentum. In the case 
of x = 1/8, two or more pre-existing haloes can merge 
with each other many times. We find that the number of 
occurrences of the merger event is nearly proportional to 
O.OIM 1 ' 25 in the case of n = and 0.02A/ in the case of 
n = —2 for larger haloes than ~ 100 cells. Then haloes with 
a larger mass generally experience many merging events. 
Moreover, the haloes with larger mass have larger angular 
momentum. Hence the contribution of the orbital angular 
momentum of a halo is greater for the halo with a larger 
total angular momentum. 

5.5 Spin parameter A 

In this subsection, we show the distributions of the dimen- 
sionless spin parameter A, which is defined as 

A " i|£|1/2 (24) 

where E is the total energy of a halo with mass M and 
angular momentum L. We estimate E by summing kinetic 
energies and gravitational energies of cells within each halo. 

This spin parameter corresponds to the ratio of the ro- 
tational kinetic energy to the gravitational energy, so an 
object with A ~ 1 is fully rotationally supported. It has 
been reported that typically A ~ 0.08 for the theoretical 
predictions and for the N-body simulations, and A ~ 0.5 for 
spirals and A ~ 0.05 for ellipticals by observations. We need 
to know dissipational gas collapse process in order to com- 
pare the spin parameters obtained theoretically with those 
obtained by observations, but it is beyond the scope in this 
paper. 

In Figjl^, we show the distributions of A in the merger 
(solid lines), homogeneous (dashed lines), and non-orbital 
cases (dotted lines) when x = 1/2. Similarly, in FigJIJ, we 
show the same lines as Figjl^ in i = 1/8. In both cases, 
the values of A at the peaks of distributions in the merger 
case are lower than those in the homogeneous case, and are 
higher than those in the non-orbital case. These tendencies 
are the same as the re sults in the previous subsections. The 
rms values of A are \/ (A 2 ) ~ 0.07 in the case of n = and 
0.25 in the case of n = —2 in the merger case. These are 
almost consistent with the previous theoretical predictions. 



Figure 10. The same as FigJ9 but for x = 1/ 



© 0000 RAS, MNRAS 000, 000-000 



10 M. Nagashima & N. Gouda 




log(L) log(L) 



1 1 1 1 1 1 1 
- (b)K-1/2 


J. ! MRG - 




r\\ MRG— 2 - 




\v\ NOR - 




\\vy NOR-2 - 


I ! /I.' / I I I I 





log(L/M 5/3 ) 



Figure 11. (a) Angular momentum distributions in the case of 
x = 1/2. The thick lines denote the case of n = 0. and the thin 
lines n = —2. The solid lines are the same as FigJq (the merger 
case), and the dashed lines denote the non-orbital case (angular 
momentum without the orbital angular momentum), (b) Mass- 
weighted angular momentum (L/M 5 / 3 ) distributions in the case 
of x = 1/2. Types of the lines are the same as (a). 



6 CONCLUSIONS AND DISCUSSION 

In this paper we have analysed the acquisition and distri- 
bution of the angular momentum of proto-objects with in- 
homogeneous density distribution by using the merging cell 
model. The assumptions in our calculations are as follows: 
spherically symmetric collapse of each block, the overlapping 
criteria (see Section 3), and the growth of angular momen- 
tum proportional to the cosmic time t. We also assume that 
the halting time of the acquisition of angular momentum is 
determined by the maximum expansion time. 

In almost all of the previous work, it has been suggested 
that proto-objects are homogeneous density ellipsoids, and 
that the maximum expansion time may be estimated from 
the smoothed density contrast. However, if we consider the 
hierarchical clustering scenarios on which most of the previ- 
ous analyses of the angular momentum are based, we should 
consider the merger effect. In the hierarchical clustering sce- 
nario, it is natural that each proto-object has a different 
halting time in the acquisition of the angular momentum. 
Hence it is very important to consider the merging histo- 
ries of dark haloes. However, we find no significant differ- 
ences between the distributions of angular momentum for 
the merger and the homogeneous cases in our analyses by 
using the merging cell model. We find the reason for this by 
dividing the angular momentum into two components, that 
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Figure 12. The same as Fig.|lT| but for x 



is, the time-independent vector and the maximum expansion 
time (see Section 5.3). The distribution of the maximum ex- 
pansion time, which is directly influenced by the merging 
history, moves to the lower value. However, the angular mo- 
mentum components that have a direction which is the same 
as that of the final angular momentum grow later than those 
with other directions. This is shown by comparing simply av- 
eraged maximum expansion time in each halo (i*' flG ) with 
the direction- weighted maximum expansion time (i m ) C osfl. 
We find that many haloes satisfy (t m ) cos g > (t^ RG ). Thus, 
as these two effects cancel each other out, the difference be- 
tween the merger case and the homogeneous case becomes 
negligible. Note that this conclusion is independent of the 
overlapping parameter x. 

In our calculations, we use the scale-free power spectra, 
P(k) oc k n , where n = and n — —2. When we normalize 
P(k) at the present epoch, the mass-scale in our investiga- 
tion corresponds to the clusters of galaxies. However, we can 
scale the normalization to a galactic scale at high redshift 
because of the scale- free spectra (see Section 5). When we 
consider the galaxy formation processes, it may be impor- 
tant to analyse the process of the acquisition of the angular 
momentum of dark haloes at high redshift corresponding 
to the galaxy formation epoch, because the spins of dark 
haloes convert to those of the baryonic component, namely, 
galaxies at high redshift. 

We also show the contributions of the orbital angular 
momentum of pre-exiting haloes when mergers occur. This 
analysis is enabled by using the semi-analytical models like 
the merging cell model. The difference between the merger 
case and the non-orbital case, which does not include the 
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Figure 13. A distributions in the case of x = 1/2. (a) n = 0. (b) 
n = —2. The solid lines are the merger case, the dashed dashed 
lines denote the homogeneous case, and the dotted lines show the 
non-orbital case. 

orbital angular momentum component strongly depends on 
the overlapping parameter x. In the case of x = 1/2, in which 
there are not too many merger events between two or more 
pre-existing haloes, the contribution of the orbital angular 
momentum to the total angular momentum is small. On the 
other hand, in the case of x = 1/8, the orbital angular mo- 
mentum plays an important role. The contribution of the 
orbital angular momentum is about one order of magnitude 
larger than the contributions of the other components, which 
are the spin angular momentum and the angular momentum 
of the accreting matter. Note that it is easy for blocks to col- 
lapse in this case of x = 1/8 and that there are many merger 
events. We cannot decide the value of the overlapping pa- 
rameter x without a detailed study of the overlapping con- 
dition, so we can only mention our conclusions qualitatively. 
However, even under this limitation, we find that the orbital 
angular momentum plays a significant role in the total an- 
gular momentum if merger events occur many times. When 
the evolution of the angular momentum of the baryonic com- 
ponent is considered, the process of the merging of galaxies 
as gaseous and stellar components will involve a release of 
angular momentum to dark matter and running away gas. 
So the distribution of the angular momentum of the gaseous 
component may follow the distribution without the orbital 
angular momentum. If elliptical galaxies are formed through 
such a merger process, and if spiral galaxies are formed with- 
out a loss of angular momentum of accreting gas, this result 
is very suggestive, because the observational difference be- 
tween the angular momentum of ellipticals and spirals is also 
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Figure 14. The same as Fig.|l3| but for x = 1/8. 



about one order of magnitude (see e.g. Fall 1983 or Fig. 9 in 
Catelan & Theuns 1996a). Such tendencies are presented in 
Fig.^. We need hydro-dynamical simulations in order to de- 
cide whether the distribution of the angular momentum of 
galaxies is followed by the distribution that includes orbital 
angular momentum or not, since our model considers only 
dark matter component. One such baryonic effects has been 
discussed as the 'angular momentum catastrophe' in Weil, 
Eke & Efstathiou (1998). 

As mentioned above, we assume the spherically sym- 
metric collapse of each block in the merging cell model. For 
this reason, we need the overlapping parameter x. So we be- 
lieve that we need a new model which includes nonspherical 
collapse in order to discuss the merger effect quantitatively. 
This study is in progress. 
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